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Abstract. We propose a phenomenological generaliza- 
tion of the models of large-scale structure formation in the 
Universe by gravitational instability in two ways: we in- 
clude pressure forces to model multi-streaming, and noise 
to model fluctuations due to neglected short-scale physical 
processes. We show that pressure gives rise to a viscous- 
like force of the same character as that one introduced 
in the "adhesion model" , while noise leads to a roughen- 
ing of the density field, yielding a scaling behavior of its 
correlations. 
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1. Introduction 

Analytical approximations for the evolution of large-scale 
structure (LSS) are based on the paradigm that small ini- 
tial perturbations grow by gravitational instability, which 
is in turn implemented in the simplest matter model, 
'dust' (Peebles 1980, Zcl'dovich & Novikov 1983, Sahni 
& Coles 1995, and ref. therein). However, this approxima- 
tion has some limitations: one has to restrict the applica- 
tion to the early stages of structure formation and when 
the effects on the evolution of physical processes different 
from gravitational instability are negligible. In this paper 
we purport to generalize this matter model in order to 
overcome some of these limitations. 

One of the problems at the later stages of LSS evolu- 
tion is the formation of multi-stream regions, i.e., regions 
where particles of dust come together with very different 
velocities. This fact manifests itself as the emergence of 
caustics in the density field, where the velocity field is 
"vertical" (i.e., where it acquires an infinite derivative at 
a point) and later on multiply valued. This problem arises 
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from insisting on following the trajectory of each parti- 
cle of dust. We therefore propose a set of hydrodynamic- 
like equations for the coarse-grained fields, which trace 
the average motion rather than that of individual parti- 
cles. A substantial ingredient of this approach is that the 
coarse-grained velocity evolves under the combined action 
of gravity and pressure-like forces due to velocity disper- 
sion (i.e., because particles of dust do not move exactly 
with the coarse-grained velocity). We cannot resort to hy- 
drodynamic considerations of local equilibrium but make 
instead use of 'equations of state' as phenomenological 
matter models without any further justification. In these 
models the pressure is assumed isotropic and may only 
depend on the (coarse-grained) density, that is, p = p{g). 
This assumption makes the problem accessible to analyti- 
cal study and helps to illuminate our major argument that 
the presence of forces which counteract the gravitational 
attraction is a basic step in understanding the dynamics 
of self-gravitating matter; a detailed study of the origin 
and properties of pressure forces is carried out elsewhere 
(Buchert & Domfnguez 1998). 

That the putative problem is far more complex than 
the dust case already becomes clear in the investigation of 
the one-dimensional Euler-Newton system with the sim- 
ple matter model p oc g (Gotz 1988). Gotz has shown that 
solutions to the one-dimensional problem can be gener- 
ated by solutions of the Sine-Gordon equation. This well- 
studied equation has a rich spectrum of solutions that in- 
cludes solitons. Gotz also pointed out that an asymptotic 
N-soliton state is generic, i.e., will be realized almost inde- 
pendently of the initial data. We see already in this com- 
paratively simple case, that we are faced with a generic 
picture which is completely different from what emerges 
in a cosmology based on dust matter: special nonlinear fea- 
tures build up structures at large times which are absent 
in the dust cosmology. This illustrates that the complex- 
ity introduced by a pressure term could bear far-reaching 
surprises. We also want to stress that the introduction of a 
pressure term is not in contradiction with our present un- 
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derstanding of LSS formation, since it partly arises from 
N-body simulations of the structure formation process 
which capture multi-streaming effects, and hence are not 
constrained by the analytical approximation to dust mat- 
ter]^ which, in its simplest realizations (e.g., Zel'dovich's 
approximation (Zel'dovich 1970)), features immediate de- 
cay of structures after their formation. 

We also consider an extension of the model to include 
stochastic effects. This "stochasticity" arises from the ef- 
fect on the dynamical evolution of physical processes oc- 
curring on time- and/or length-scales much smaller than 
those directly associated with LSS formation, thus allow- 
ing to model them by means of a stochastic source (a 
noise). Possible sources are deviations from the mean field 
approximation, fluctuations inherent to the hydrodynamic 
(i.e., coarse-grained) description, and non-gravitational 
processes in baryonic matter. We shall use the simplest 
model of a Gaussian-distributed forcing on the coarse- 
grained evolution. 

As with pressure-like forces, we just want to stress that 
a noisy forcing could be relevant to LSS formation, but a 
detailed consideration of its origin and properties is be- 
yond the scope of the present paper. In fact, application 
of the Renormalization Group shows that noise is rele- 
vant in non-exceptional conditions (Barbero et al. 1997; 
Dommguez et al. 1999), implying that even if it is very 
weak (apparently negligible), its effects are amplified and 
can have a non-negligible effect on the dynamical evolu- 
tion. 

This paper is structured as follows: in Sect, ^we begin 
by presenting the basic system of equations in the Newto- 
nian regime, and then proceed to a discussion of restrictive 
assumptions for the weakly nonlinear regime. In Sect. I we 
discuss the role of the pressure-like force for some partic- 
ular choices of the equation of state p = p{g) and the 
connection with Burgers' equation. In Sect, ^we consider 
the role of noise and provide a detailed description of the 
relationship between the cosmological equations and the 
Kardar-Parisi-Zhang (KPZ) equation. In Sect. |^ we study 
the linear regime in the presence of pressure and noise. We 
finally conclude in Sect. Some technicalities have been 
left for two appendices, one devoted to the exploration of 
the validity of what we call in Sect. ^ the "adiabatic ap- 
proximation" , and another to a more detailed discussion 
of the linear regime. 



2. Basic Equations and Restrictive Assumptions 

We are interested in discussing LSS formation in the non- 
relativistic regime and therefore consider the Newtonian 
cosmological equations for a self-gravitating fluid in a 



standard Friedmann-Lemaitre (FL) cosmological back- 
ground dominated by non-relativistic matter (Peebles 
1980). The cosmological background is characterized by 
the cosmic expansion factor a{t) and the homogeneous 
background matter density Qb{t), which obey (Bubble's 
function is defined as iJ = a /a) 



, SttG K a 



Qb = Qoa 



(1) 



where the constant K determines the sign of the spatial 
curvature (treated as an integration constant in the New- 
tonian framework considered throughout this paper), Qq 
is the background density at some time to when a (to) = 1, 
and A is the cosmological constant. Without loss of gen- 
erality, one can choose to to correspond to the present 
epoch. 

It is convenient to work in comoving coordinates 
X = a^^r, where r are the standard non-rotating Eule- 
rian coordinates. The fundamental fields will be as fol- 
lows: the density g (or equivalently, the density contrast 
S := [qI Qb) — 1 ), the peculiar-velocity u = VLphys — Hr, 
where Uphys is the physical velocity and Hr is the Hub- 
ble flow, and the gravitational peculiar-acceleration g — 
Sphys + 5 (47rGpb - A) r, where gphys is the physical grav- 
itational acceleration and — ^ {AnGgb — A) r is the Newto- 
nian counterpart of the gravitational acceleration oppos- 
ing to the background expansion. We subject g, u and g to 
periodic boundary conditions on some large scale to assure 
uniqueness of the cosmological solutions, in which case gb 
is equal to the spatially averaged density (see Ehlers & 
Buchert 1997). 

The fields p, u obey a set of hydrodynamic equations 
expressing the conservation of mass and momentum in 
an expanding background. These equations are quite 
similar to those of the standard dust model (Peebles 
1980), except for two forcing terms in Euler's equation 
that model the multi-streaming and stochastic effects 
discussed in the Introduction 0: 



• Continuity equation: 
dg 



, 3Hg+-\/-igu)^0; 
at a 



(2) 



• Euler's equation: 
du 1 



1 



^ We here imply 'single-streamed dust' as opposed to models 
which also cover multi-dust regions. 



, (u-V)u + iJu = g Vp + s; (3) 

at a ag 

• Newtonian field equations: 

V • g = -AirGaig ^ gb) , V x g = . (4) 

We emphasize that the integral curves of the peculiar- 
velocity field u are not associated to trajectories of in- 
dividual particles, rather g, u and g are considered as 

^ From now on, time derivatives are taken at constant x and 
gradients refer to comoving coordinates. 
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coarse-grained fields. This coarse-graining is the origin of 
the two new terms on the right-hand-side of Euler's equa- 
tion (|^) . One of these new terms is the pressure force Vp, 
which accounts for the isotropic part of the multi-stream 
force, and therefore models velocity dispersion (that is, the 
fact that in any infinitesimal cell there are particles with 
different velocities). Because of this, the integral curves 
of u represent trajectories of the mean (possibly multi- 
streamed) flow after averaging over velocity space. This 
pressure term is not related to thermal pressure, which 
can be indeed neglected on the scales we are interested in. 
It is a model of the velocity dispersion generated by grav- 
itational instability (see Buchert & Dommguez 1998 and, 
for a recent generalization to general relativity, Maartens 
et al. 1999). 

The other new term is the stochastic force represented 
by the noise s; it accounts for processes hidden by the 
coarse-grained description of the fluid and whose typical 
time- and length-scales are much shorter than those ex- 
plicitly considered for LSS formation. We have resorted 
to modelling these processes as a stochastic forcing and 
include: (a) the effects of small-scale degrees of freedom 
whose physics is also governed by non-gravitational pro- 
cesses, (b) deviations from mean field behavior, manifested 
as random forces acting on the particles of the gravita- 
tional gas as a consequence of independent impulses of 
random size and amplitude arising from "sling-like" pro- 
cesses in encounters (Kandrup 1980), (c) deviations of the 
density and velocity fields from the values prescribed by 
the deterministic version of Eqs. (|[^) due to the graini- 
ness of the underlying physical system of particles (Lifshitz 
& Pitaevskn 1980). 

To close the system of Eqs. (|[^) a relation is needed 
between the dynamical pressure p and the two indepen- 
dent fields g and u, as well as a specification of the statis- 
tical properties of the stochastic force s. As for the former, 
we assume the local relationship p = p{g). There does not 
seem to be any a priori reason for this "slaving" of pressure 
to density (one cannot resort to the hypothesis of local 
equilibrium as is done for the thermodynamical pressure 
in fluids), so the success of this phenomenological assump- 
tion must be judged according to the conclusions following 
from it. In fact, a detailed study of the origin of pressure 
forces in Eq. (||) provides p oc g^^^ under the assumption 
of small velocity dispersion (Buchert & Domi'nguez 1998) 
and therefore p = p{g) is the most straightforward phe- 
nomenological generalization. We also require p'{g) > 0, 
that is, pressure opposes gravitational collapse. 

As for the noise, we make the assumption of Gaussian 
distributed noise. Since noise is due to the short-scale de- 
grees of freedom, this assumption could be justified by the 
central limit theorem. As is well known, Gaussian noise 
can be characterized by just two moments: its mean, which 
we require to vanish, (s) = (since any systematic forc- 
ing should be made explicit in Eq. (§)), and its two-point 



correlations 



(s,(x,i)s,(x',t')) = 2Aj(x,x',t,f') 



(5) 



where (x, x', t, t') is the covariance matrix with mixed 
discrete and continuous indices (Gardiner 1994; Van Kam- 
pen 1992). We ignore the possibility that Dij depends on 
g or 11 (quenched noise), since this renders the analysis 
of Eqs. (|^^ too difficult. As with pressure, this must be 
seen as a phenomenological assumption, whose merit will 
be judged from the final results. Notice however that we 
consider the possibility of colored Gaussian noise, i.e., the 
noise s is correlated over space and time, as expressed 
by the (in general) non-trivial dependence of D on posi- 
tion and time in Eq. (|^) . Colored noise represents a richer 
texture of physical effects than white noise and the partic- 
ular case of power-law correlated noise is still amenable to 
analytical study by means of the Renormalization Group 
(Dominguez et al. 1999). We later restrict the generality 
further by choosing a curl-free noise and characterize it 
statistically by a single function D instead of Dij. 

The analytical study of the system of Eqs. (|^^ is 
very difficult in general. Inspired by the deterministic dust 
case (i.e., p = 0, s = 0), we simplify the problem and 
put forward the assumption of parallelism: we impose the 
condition that the peculiar-velocity is a potential field and 
remains parallel to the gravitational peculiar-acceleration 
field: 



g = F{t)n 



(6) 



where F{t) > is a proportionality coefficient which fol- 
lows from the deterministic linear theory for dust (see 
Sect, ^and in particular the discussion after Eqs. (!65[-^66|) 
in Appendix B). This assumption implicitly requires that 
both pressure and the intensity of the noise are small when 
compared with the dominant self-gravity in Eq. (^), which 
restricts our considerations to spatial scale regimes close 
to the validity limit of the dust model; noise and pressure 
will typically dominate on scales small compared to this 
limit. The assumption of parallelism underlies the well- 
known "Zel'dovich approximation" (Zel'dovich 1970), and 
is well-justified for deterministic dust models in the lin- 
ear as well as in the weakly nonlinear regimes (see Bild- 
hauer & Buchert 1991; Kofman 1991; Buchert 1992; Hui 
& Bertschinger 1996; Susperregi & Buchert 1997). This 
(indeed oversimplifying assumption) will be very useful 
to analytically access the problem, and to define "local" 
approximations. 

Also, we study that assumption first because popular 
models like the "adhesion approximation" (Gurbatov et 
al. 1989) can be derived on the basis of this assumption. As 
we shall see, this assumption is consistent with the picture 
that, coming from large scales, the mean motion is ruled 
by the dust model, but incorporation of the effects of pres- 
sure and noise adds several interesting aspects to it. At the 
smaller scales, the "parallelism assumption" requires that 
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the "backreaction" of these effects on the trajectories of 
the mean dust flow be neglected. In the dust case it can be 
shown that this assumption admits a class of 3D solutions 
(Buchert 1989); this class is highly restrictive, but using 
it as the basis for approximation schemes turns out to be 
very successful (e.g. Buchert 1996 and ref. therein) - a 
subcase of these schemes is known as Zel'dovich's approx- 
imation. However, in the present case we cannot disprove 
that the assumption (H) is too restrictive to allow for the 
existence of a class of exact solutions. We will learn more 
about the justification of this assumption in Sect. ^. 

Under the assumption (^, Eq. (||) reduces to 



du 1, 



V)u- 



ag 



(7) 



with constraints following from the field equations (|j) : the 
velocity field is irrotational, V x u = 0, and its divergence 
satisfies 



(8) 



The pressure force term Vp can be computed from this 
last expression, which gives 



Vg 



iirGa 



and Euler's equation may now be finally written as: 



1 , 

— + - u- V u- 

at a 



{H - F)vL 



(9) 



(10) 



where we have defined a coefficient v which behaves as 
a kinematic viscosity and which we call "gravitational 
multi-stream" (GM) coefficient (because it has its ori- 
gin in the interplay between self-gravitation and multi- 
streaming flow); it is given by 



AnGaityg 



> 



(11) 



and depends on the density, and explicitly on time through 
F{t) and a(t). As we see from Eqs. ( p^ ) and (pd|), the 
pressure behaves in Euler's equation effectively as a vis- 
cous force that prevents caustics in the velocity field, a 
behavior that essentially requires the participation of self- 
gravity. The interplay between self-gravity and pressure 
leads to the stabilization of the structures, and holds them 
together. We also want to stress that the GM coefficient 
does not generate dissipation, since our starting set of Eqs. 

lacks any sort of dissipative term (see Buchert & 
Domfnguez 1998 for further details about this apparently 
paradoxical point). 

The last step is the elimination of the density g in the 
expression of the GM coefficient in favor of V • u by means 
of the constraint (||), thus reducing Eq. ( [l0| ) to a closed 
equation for the velocity field u. This equation will be 
explored in the next sections. 




Particle displacement 

Fig. 1. The coarse-graining Idea is exemplified in a 
schematic way based on a 2D tree-code simulation 
(Buchert 1996): the particles are attracted towards the 
pancake, a multi-stream region containing many streams 
(velocities) at a given Eulerian position. The coarsening 
cell follows initially a trajectory similar to that of individ- 
ual particles. But as the cell moves into the pancake, the 
kinetic energy of the bulk motion is gradually transformed 
into internal kinetic energy, and its trajectory becomes 
qualitatively different from that of a particle. 



3. The Role of Pressure: Burgers like Equations 

In this section we consider the role of pressure in the dy- 
namical evolution and thus for the time being we drop the 
noise term, s(x, t). But before going into the details, it is 
worthwhile to get an intuitive picture of how the pres- 
sure term affects the dynamical evolution. As remarked 
above, pressure has in Eq. (|l^) a stabilizing effect oppos- 
ing collapse. But this effect holds even in regimes where 
the assumptions leading to Eq. (|l^) (especially the paral- 
lelism assumption) break down: Fig. yj shows the evolution 
of coarsening cells in a simulation and how velocity dis- 
persion influences their motion. 

Different relationships between pressure and density 
yield different ^^-dependences of the GM coefficient. In this 
section we will concentrate on polytropic models p = K.g'^ , 
where n is constant and 7 is the polytropic index. The 
prefactor k should be small enough so that the paral- 
lelism assumption (^ is approximately satisfied on the 
length scales we are interested in, and it should also sat- 
isfy K7 > in order to fulfill the condition p'{g) > 0. The 
choice of polytropic models is interesting because we are 
able to recover well-known approximations in cosmology 
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(Zel'dovich's approximation, the "sticky particle model" 
and the "adhesion model") as well as other interesting 
models, and we can understand them within a single equa- 
tion. With this form of p{q), Eq. (nn) becomes 



^ + i(u.V)u+(i7-F)u: 



rz-^V^u. (12) 



We employ now Eq. (||) to eliminate g in place of V • u, 
and since V^u = V(V • u) (because V x u = 0) one gets 



O 1 

-^ + -(u-V)u+(H - F)u = -i>V[l - /3V • uP-^ 
at a 



(13) 



where v = ^k,qI ^ / {l 



l)a and [3 = F/AnGagi, are func- 
tions which depend only on time. This equation can be 
further simplified as follows: making explicit use of b(t), 
the growing mode of the density field in the linear regime 
of the deterministic dust case, one has F = AirGgbb/b (see 
Eqs. (^6^.66) in Appendix B and discussion thereafter). 
Hence, defining a new velocity field v = u/ab and employ- 
ing b as the new time variable, Eq. (O) can be written as: 



^ + (v V)v = -/xV[l - 6V • v]^-i 



(14) 



When 7 = 2, this equation becomes the 3D Burgers' equa- 
tion (Burgers 1974) (except for the time-dependence of 
the GM coefiicient /i). In the cases 7 ^ 2, we are dealing 
with generalizations of Burgers' equation. We now discuss 
several cases for the choice of 7 ^ 

Zel 'dovich 's approximation 
If one chooses 7 = 0, then /i = and Eq. (|l^) reduces 



to 



dv 

- + (v.V)v = 



(15) 



which corresponds to the well-known Zel'dovich's approx- 
imation; that approximation is equivalent to a subclass of 
Lagrangian first-order solutions, the subclass being just 
defined by our assumption (^ (see Buchert 1989, 1992). 
The dynamical evolution governed by this equation gener- 
ically leads to the formation of singularities where the ve- 
locity field is multi-valued (Arnol'd et al. 1982). 



^ Eq. (^3]) is valid for 7 7^ 1. If 7 = 1, one has a logarithmic 
function instead of a power on the right-hand-side of Eq. (|l^), 
see Eq. 

* The names given below to some of the models stem from 
thermodynamical notions; we kept those names but do not 
imply that we describe situations in thermodynamical equilib- 



Adhesion approximation 

With the choice 7 = 2, the GM coefficient is inde- 
pendent of the density, and Eq. ( p^ ) becomes Burgers' 
equation with a time-dependent "viscosity" : 



+ (v V)v = fibV^v 

00 



(16) 



This equation is formally equivalent to the so-called "ad- 
hesion approximation" (Gurbatov et al. 1989; Weinberg & 
Gunn 1989), aside from the time-dependence of the GM 
coefficient Q 

Letting k — > (which implies /i — > 0) in this model, 
we get "maximal adhesion" and the large-scale structure 
is built from a skeleton of the "honeycomb-type" : the 
evolution is governed by Zel'dovich's approximation, Eq. 
(|l5|), everywhere except at the singularities, which become 
shock fronts. This limiting model is known as "sticky par- 
ticle model", for which geometrical contruction methods 
have been developed (Pogosyan 1989; Kofman et al. 1992; 
for review see Sahni & Coles 1995). In our derivation this 
limit implies p — > 0, i.e., vanishing velocity dispersion. It 
must be noticed, however, that this limit is singular: if 
p = strictly, then we would recover Zel'dovich's approx- 
imation, which is qualitatively different from the "sticky 
particle model" . Therefore, a non- vanishing, though very 
small, velocity dispersion is required to recover typical 
properties of the "adhesion model" . 



Isothermal model 

This model corresponds to 7 = 1 and is mostly stud- 
ied in connection with the linear theory of gravitational 
instability (since linearization of the pressure term in Eq. 
(0) amounts to choosing 7=1). Eq. (|2|) yields: 



dv 
'db 



+ (v ■ V)v 



a^b^ 



Vln(l - 6V • v) 



(17) 



It is worth mentioning that Eq. (|T^) particularized to 7 = 
1 resembles the Navier-Stokes equation with respect to 
the density dependence of the GM coefficient. 

As pointed out in the introduction, Gotz (1988) has 
shown that the asymptotic state of this model in one- 
dimensional space and neglecting the expanding back- 
ground is a collection of solitons. On the other hand, par- 
allelism is not an approximation but an exact relationship 
in one dimension (although the proportionality coefficient 
will be in general a function of time and position). Hence, 
one can expect that Eq. (17) leads to a similar picture to 
that obtained from the adhesion model, which is in turn 



^ In an Einstein-de Sitter background cosmology bfi ~ b~'^ 
and thus the GM coefficient in (llffl approaches zero as time 
goes by. Remember also that in the "adhesion approximation" 
the constant coefficient was introduced ad hoc to phenomeno- 
logically model gravitational "sticking of particles" as observed 
in N-body simulations. 
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consistent with that offered by Gotz: the shock fronts, sta- 
bihzed by velocity dispersion, will play the role of solitons. 



Adiabatic model 
For this model 7 = 5/3 and Eq. becomes 

1^ + (v • V)v = - bV • v]2/3 . (18) 

The interesting point about this choice is that there are 
physical reasons to prefer this model over any other one in 
the early stages of LSS formation: Buchert & Domi'nguez 
(1998) show how this form of p{q) can be derived from 
dynamical considerations. 



Cosmogenetic model 

By choosing 7 = 4/3 we recover the "cosmogenetic" 
model (Chandrasekhar 1967): 



9v 
'db 



(v-V)v = -^V[l-6V-v]i/3 



(19) 



This model is of interest because it is the only polytropic 
model compatible with "comoving hydrostatic equilib- 
rium", i.e., with the solutions of Eqs. (|2[^) corresponding 
to u = 0, so that the temporal evolution of the gas simply 
follows the expansion: Eq. (||) yields g{x, t) = a~'^ it) gQ{~K) , 
Eq. (^) then becomes g(x, t) = a~^(i)go(x). The "hydro- 
static equilibrium" condition follows from Eq. (^ (drop- 
ping the noise): Vp = agg. Combining this condition with 
the polytropic relationship yields 



-fna '-^'^qI VfiQ 



£'ogo 



(20) 



which can be satisfied only if 7 = 4/3. (Obviously, the 
state of "comoving hydrostatic equilibrium" itself is in- 
compatible with the parallelism assumption). 



To conclude this section, one can conjecture that 
Burgers-like equations ([l^ ) lead quite generically in the 
limit K — > to the same picture: the dynamical evolu- 
tion would be governed almost everywhere by Zel'dovich's 
approximation except at the caustics, where the right- 
hand-side of Eq. (|lj) would dominate the evolution and 
a shock structure would be formed (known as a "pancake" 
in the cosmological literature) . Only the details of the den- 
sity and velocity profiles in the neighborhood of pancakes 
would depend on 7. It is an open question whether and 
in which cases velocity dispersion and self-gravity could 
balance to form stable soliton-like configurations. 



4. The Role of Noise: Emergence of a KPZ like 
Equation 

In this section we consider the dynamical effect of the 
noise. We come back to Eq. (0) and introduce a veloc- 
ity potential ip{x.,t) by u = — VV' (since the parallelism 
assumption (0) implies that u is irrotational) . Then 



'4 - 

ot 2a 



[H - F)ij 
1. 



'dym 

Q, y 



(21) 



We now make the assumption that the noise s is a poten- 
tial forcing, i.e., that there exists a (stochastic) potential 77 
such that s = — V77. This stochastic potential will be cho- 
sen Gaussian-distributed with zero mean, irf] = 0, and 
correlations 



(7](x,i)77(x',i')) = 2Z?(x,x',i,i') 



(22) 



The assumption of potential noise can be motivated by the 
fact that in the linear regime only the potential component 
of the noise contributes to the growth of perturbations by 
gravitational instability (see Sect. ^ and Eq. (.58) in Ap- 
pendix B). If the noise had a non-potential component, it 
would generate vorticity in the velocity field u and would 
invalidate the parallelism assumption. 

After inserting this stochastic force into Eq. (|2^) , inte- 
grating, dropping an irrelevant additive function of time 
and changing from gi to (5 in the integral, one obtains: 



1 



9-0 



dt 2a 



{Vipf + {H- F)ij 



dy 



p'{gb{l + y)) 
1 + y 



+ V- (23) 



On the other hand, the constraint (H) can be converted 
into an equation relating the velocity potential ^ and the 
density contrast 6, that is 



5(x,i) = 



Fit) 



ATrGa{t)gbit) 



V20(x,t) . 



(24) 



Inserting this equation into Eq. ( |23|) above then yields 
a non-linear partial differential equation for the velocity 
potential (after specifying the matter model p{g)) given 
by: 

f -^(V0)^ + (i/-^)^ = 

1 f^-'\^PMl±y)l^,. (25) 



l + y 

There are two sources of non-linearities: the convective 
(also known in the literature as advective) term (VV')^ and 
the integral arising from the pressure force. To proceed 
further we need an explicit expression for p{g) to evaluate 
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the integral. For a polytropic model p = Kg'^ , one recovers 
Eq. but expressed in terms of the potential ip (plus 
the noise term). As discussed at the end of the previous 
section, one may expect that the gross features of the LSS 
emerging from these models are quite insensitive to the 
particular equation of state p = p{g) in the limit of small 
pressure and that only the fine details depend on it. We 
therefore simplify this equation further by expanding the 
integral in a Taylor series and keeping only the lowest 
order term in 6. The resulting equation is exact for the 
particular choice p = kq^: 



(«-F)* = ^5g!^vV + ,. (26) 



As remarked in the previous section, this choice leads to 
the "adhesion approximation" , which is in fact postulated 
on a similar reasoning: it is the simplest way of mod- 
elling "sticking particles" (Gurbatov et al. 1989). Hence, 
the "adhesion approximation" can be viewed as the first 
term in a Taylor series that approximates the integral ex- 
pression. 

Eq. (|2^) is the simplest equation for ip that we can 
write, still containing the main ingredients that enter into 
the physics of the self-gravitating gas we are describing, 
although as we will see, it already entails a considerable 
degree of complexity. 



Eq. describes what is known in the condensed 

matter literature as interface growth phenomena, but with 
two added ingredients which are substantive to cosmology: 
the presence of time-dependent coefficients due to cosmo- 
logical expansion, and the presence of a term proportional 
to V', which in the context of condensed matter physics is 
interpreted as a finite, albeit time-dependent, correlation 
length. This equation (Berera & Fang 1994; Barbero et al. 
1997) is then a generalization to cosmological settings of 
the KPZ equation (Kardar et al. 1986, Barabasi & Stanley 
1995) for surface growth. In principle, one could use tech- 
niques similar to those used there to study this equation, 
but it turns out that complications arise as a consequence 
of the time dependence of the coefficients. Because of this, 
the first thing one thinks of is to perform changes of vari- 
ables that will bring the equation into an equation with 
constant coefficients from which one can later on proceed 
with the analysis. With this in mind we rewrite Eq. (26) 



T 



^ + 77, (28) 



where we have defined the following dimensionless func- 
tions of time: 



1 F{t)p'(Mt)) 



V AT:Ga?{t)Qb{t) ' 



(29) 



4.1. Discussion of Eq. 

This stochastic partial differential equation is first order in 
time and second order in position; it is also non-linear in 
^1) and the coefficients of the different terms are functions 
of time. The physical meaning of each term is transparent: 

(a) The term proportional to ijj encompasses the compe- 
tition between damping of perturbations due to the cos- 
mological expansion, (Hip), and enhancing of perturba- 
tions due to gravitational collapse, {—F^j). This term in- 
troduces a time-dependent time scale \H — F\~^ ^ which is 
the time scale for the damping (or enhancing) of pertur- 
bations in regimes when the nonlinearity is negligible. 

(b) The term proportional to the Laplacian describes the 
damping of perturbations due to velocity dispersion. This 
term defines a time-dependent length scale, Jeans' length. 



p'iQb) 



47rGa2 



Qb 



1/2 



(27) 



which is discussed in Appendix B, after Eq. (.63). 

(c) The non-linear term is the convective term in Euler's 
equation (^. The effect of this term is to broaden the 
peaks in the field ip. (An intuitive picture of how a term 
like (V'0)^ behaves may be found in Barabasi & Stanley 
1995, Fig. 6.2). 

(d) The noise term incorporates the effects of fluctuations 
due to different sources and leads to a roughening of the 
field ip in space and in time as it evolves. 



hit) 



\a{t) ' 



h{t) = [F{t)-H{t)]T . 



(30) 



(31) 



The dimensional parameters v, A, T are introduced as 
bookkeeping quantities to carry the dimensions. 

Defining a new time coordinate r, a new velocity po- 
tential ^'(x, r) and a new noise ^(x, r) via 



T{t) - / dyhiy) , * = , ^ = 4'7 , 

Jta /l /l 

allows one to recast this equation into: 



(32) 



^^yVH + hsj-^f^r{T)m+^, (33) 

OT Z 



where 
r(T) 



dfMr)) 



dh{t{T)) 



h{t{r)) dr 



/2(i(r)) dr 

1 /3(i(r)) 



Tfi{t{r)) 



(34) 



The above equation is a "massive" KPZ equation, but 
with the peculiarity that the coefficient of the term pro- 
portional to ip (the "mass" term) depends on time: one 
has a standard KPZ equation if r(r) = 0, time-dependent 
damping of the surface growth if r(T) > 0, or explosive 
unstable behavior for r(r) < 0. 
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Due to the noise term in Eq. (p3|), the field will de- 
velop dynamical correlations in addition to those due to 
the initial conditions. If the nonlinearity were neglected in 



Eq. (33), the effect of noise would simply be to superim- 
pose Gaussian fluctuations on the deterministic evolution. 
But the nonlinearity couples different length scales, so that 
the evolution at any given scale also receives a contribu- 
tion from fluctuations on other scales. Eq. (^) implicitly 
involves a coarse-grained description and thus a smooth- 
ing length scale, so that the nonlinear coupling promotes 
the phenomenological coefficients i/. A, r (and the noise 
itself) into scale-dependent quantities. 

The correlations induced by noise can be computed by 
means of the Renormalization Group (Gell-Mann & Low 
1954, Binney et al. 1993, Weinberg 1996), via the com- 
putation of the scale-dependence of the coefficients. The 
fact that r(r) depends on time again complicates the ap- 
plication of the Renormalization Group. However, if r(r) 
happens to be independent of t or if an adiabatic ap- 
proximation is justified (i.e., if r^r) varies very slowly on 
the time-scales associated with the dynamical evolution 
prescribed by (|3^)), then one can to some degree of ap- 
proximation neglect the r-dependence of r(r) when apply- 
ing the Renormalization Group. In Appendix A a detailed 
study of this question is carried out. 

Under the adiabatic assumption, application of the 
Renormalization Group is straightforward albeit analyt- 
ically cumbersome, and has been carried out elsewhere 
(Barbero et al. 1997; Domfnguez et al. 1999). We here con- 
tent ourselves by quoting the conclusions. The KPZ equa- 
tion exhibits self-affine correlations in the large-distance, 
long-time regime, and the "massive" KPZ equation, Eq. 
(p3|), can also exhibit this behavior. In particular, the 
equal-time two-point density correlation obeys in such 



case the scaling (5(x, t)5(x', t)) 



X 



'|2x-4 



where x 



is the "roughness" exponent, which can be computed by 
means of the Renormalization Group. This result shows 
that noise may be relevant in that it can induce the gen- 
eration of self-affine correlations in a self-gravitating col- 
lisionless gas, even if the initial conditions are not self- 
affine. 



5. The Linear Theory Revisited 

The strongest assumption that we have made has been the 
parallelism hypothesis (||), motivated by the well-studied 
dust case, where it is justified in the linear and weakly 
nonlinear regimes. In this section we will study the lin- 
earized version of the system of Eqs. (^^ to find out how 
well this assumption is justified in the presence of pressure 
and noise. 

The linearized Newtonian cosmological equations 
around the unperturbed FL background (i.e., ^ = 0, 
u = 0, g = 0) read: 



• Continuity equation: 

dt a 

• Euler's equation: 



V • u = ; 



^ + Hu = g Vd + s 

ot a 



Newtonian field equations: 



V • g = -AnGagbS , V x g = 



(35) 



(36) 



(37) 



Since the noise s(x,i) has vanishing mean and the 
equations are linear, the general solution can be split into 
the sum of the averaged fields {S), (u) and (g), obeying 
the deterministic version of the Eqs. (35 -^7|) (i.e., setting 
s = 0), plus a fiuctuating part which is a linear functional 
of the noise s. The dynamical evolution thus follows a de- 
terministic trajectory with Gaussian fluctuations super- 
imposed on it, but these corrections will be small if the 
noise intensity is small. The detailed study of this sys- 
tem of equations is performed in Appendix B. The main 
conclusions for our purposes are the following: 

(i) The average vorticity of the peculiar-velocity is 
damped by the background expansion, even if p 7^ 0, 
and thus (u) becomes asymptotically a potential field as 
t +00. 

(ii) The parallelism assumption (^ gets modified due 
to pressure. When the lowest order correction is included, 
we find that f\ 



(g)«Fo(t)(u) 



p'{g^)Fl{t) 



V^(u) 



(38) 



This equation shows that corrections to parallelism in the 
average motion are negligible in the limit of small pressure 
far from caustics (i.e., where (u) is a smooth field and 
V^(u) does not diverge). Pressure becomes important in 
pancakes, where parallelism no longer holds and also the 
assumptions of small and isotropic pressure break down. 



6. Conclusions 

In this paper we have put forward the set of Eqs. 
as a new, phenomenological approach to the problem of 
LSS formation. The difference with regard to similar ap- 
proaches lies in the interpretation of these equations as 
describing the dynamical evolution of the coarse-grained 
fields S, u and g, which is the origin of the pressure-like 
force and of the fluctuations (noise) . 

We have shown that under the parallelism assumption 
(^, the pressure force gives rise to a viscous-like force of 
the same character as that of the "adhesion model" , which 
is a successful model on large scales. As is also known in 

® The notation and the assumptions leading to this expres- 
sion are explained in Appendix B. 
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the context of the "adhesion model" , the limit of vanishing 
pressure is singular: the models with no pressure {p — 
exactly) are qualitatively different from the models with 
p 0. Therefore, even if pressure is very small (seemingly 
negligible), one should not set p = 0. 

We have also considered the effects of fluctuations. Un- 
der the parallelism assumption, the problem can be cast 
into a "massive" KPZ model of surface growth phenom- 
ena, which can exhibit self-afBne correlations because the 
noise can be relevant, even if it is vanishingly small. Hence, 
the limit of vanishing noise is also singular and the same 
warning as with pressure applies. 

Finally, we have explored the plausibility of the paral- 
lelism assumption (|^) by studying the linearized equations 
and reached the conclusion that it is justified far from 
caustics and close to the limit of vanishing pressure and 
noise. In fact, given the number of assumptions we have 
made, there is plenty of room to generalize the results 
presented here: relaxing the parallelism assumption, gen- 
eralizing the equation of state p = p{g), modelling multi- 
streaming by an anisotropic stress-tensor 11^ in Eq. (^ 
rather than by a pressure force (see, e.g., Maartens et al. 
1999), as well as performing a full Renormalization Group 
analysis with time-dependent coefficients. 

Furthermore, we have given the formal hydrodynam- 
ical basis on which to anchor future studies leading to 
an understanding of various scaling relations found in the 
context of LSS. 
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Appendix A: Adiabatic approximation 

In this appendix we explore under what conditions the 
coefficient r(T) in Eq. ( p3| ) can be considered almost time- 
independent. Our starting point is Eq. (p4): 



r(r) 



1 



dfMr)) 



1 



df2{t{T)) 



fMr)) dr 



/2(t(r)) dr 

fMr)) 



The function F in the definitions (|2' 
of Appendix B. Therefore: 



Tfi{t{T)) 
ih obeys Eq. 



.39) 



65D 



}_dh 
fi dt 



F 



p'iQb) 



(.40) 



Taking these results into Eq. ( |.39| ) yields 
ATiGgbO^ 



Fp'iQb) 



l!^-,2(i/-F)-3i/,.^ 

F p'iQb) 



(.42) 



On the other hand, one can identify two time scales in Eq. 
331): a time scale due to expansion 



(.43) 



and an intrinsic time scale Tint ■= kl^^- If relative 
variation of the coefficient r(T) on these time scales is 
small, then an adiabatic approximation is justified and r 
can be considered to be a slowly varying function of time. 
To quantify this statement we define adiabatic indices as 



FcKp • Fqxp 



1 dr 




1 


dr 


r dr 




Hr 


'dt 


1 dr 




1 


dr 


r dr 




fir' 


dt 



(.44) 



(.45) 



exp 5 Fi^i <^ 



The adiabatic approximation then requires 
1. We now discuss the various background cosmologies 

A.l. Einstein~de Sitter background (K 
We first study the simplest case, for which 

2/3 



a{t) 



, i^W = |, ..(^) ^ (.46) 



From Eq. ( |.68| ) of Appendix B we find that F = t ^. 
Taking these expressions into Eq. (.42) then yields 



We are interested in polytropic models p = kq'^ , for which 
the previous expression becomes 



r{t) 



4(6^G)^-i(l-7) ^|(3. 



3K7iQ^'^ 



-4) 



(.48) 



We find some interesting conclusions: for the isothermal 
model, J — 1, we get r — exactly and Eq. (|3|) reduces 
to a bona fide KPZ equation. For the cosmogenetic model 
7 — 4/3, we get r{t) = tq < independent of time. For 
other matter models we have from Eq. (1.4f; 



Eexp = |37 - 4| , Eint = 



|37-4| 
3|7-1| 



(.49) 



/2 dt 



Then the adiabatic approximation will be justified when 
(41) close to the cosmogenetic value of 4/3, in which case 

r < 0. 
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Fig. .1. Kr versus time in an undercritical FL background 
{flo < 1). Kr is measured in units of 10~'^(1— fio)^^ Gyr~^, 
and time is measured in units of r2o(l — fla)^'^/^h^^ Gyr. 

A. 2. Background cosmologies with (K 7^ 0, A = 

In this case we can also obtain explicit expressions for 
the coefficient r and the indices E^xp, Eint analytically, 
but the algebra is rather involved and the final expres- 
sions are best studied numerically. We therefore skip the 
algebra and go right to the final results, particularized to 
the isothermal model p = Kg, because this leads to the 
simplest expression in Eq. ( |.42| ). 

In an undercritical universe {K < 0) the coefficient r 
is a positive, monotonically decreasing function of time. 
In Fig. [1] we plot Kr versus t Q The adiabatic indices are 
monotonically increasing functions of time but they are 
bounded: 



32 

< Eexp < 1 , 75 — 



< 



1 



(.50) 



In the range 0.2 < ilo < 1 and at epochs t < 10h~^ Gyr 
one gets < E^^p ^ 0.35. Although the adiabatic indices 
are not much smaller than one, one can still work under 
the adiabatic assumption to a first approximation. 

In an overcritical universe {K > 0) the coefficient r 
is a negative, decreasing function of time. In Fig. we 
plot Kr versus t. The index Ei„t decreases in time and is 
bounded: 



27r2(97r2 - 64) 
(37r2 - 64)2 



< Eint < — 



32 
75 



(.51) 



The index Eg^p increases with time and is unbounded: it 
diverges at the epoch of maximum expansion because then 
the expansion rate H vanishes. This is, however, physically 



We have introduced Hubble's constant at the present 
epoch. Ho = O.lhGyr^^ , and the density parameter, Qo — 
8nGgo/3Hl 



Fig. .2. Kr versus time in an overcritical FL background 
(f2o > 1) up to the epoch of maximum expansion. Kr is 
measured in units of 10~^{^lo — l)h'^ Gyr^^, and time is 
measured in units of r2o(f^o ~ l)~^/2/j-i Gyr. 

uninteresting. At epochs t < lOh^^ Gyr and with the con- 
servative bound rig ^ 2, one finds Eint ^ 0.5 and the same 
remarks apply as in the undercritical universe. 

Therefore we conclude that within the physically in- 
teresting range of values of the parameters, an adiabatic 
approximation for the coefficient r may be justified. 



Appendix B: Linear regime 

In this appendix we study the linear regime when pres- 
sure is taken into account. The linearized system of Eqs. 
( ^5[]37| ) can be more easily studied in Fourier space. We 
therefore define the Fourier transformation (assuming a 
box of volume V): 



— ik-: 



d^x ,5(x, Oe**" " , (.52) 



and similar definitions for u(k,t), g(k, i) and s(k, i). In 
Fourier space the linearized Eqs. (pqJ3^) read: 



d6 
dt 



-ku = 



+ Hu = g + ^ ok -f s 

ot a 



ik - g — inGagbS , k x g = 



(.53) 

(.54) 
(.55) 



We decompose the velocity field in transversal u± (sat- 
isfying k • ux = 0) and longitudinal components U|| (sat- 
isfying k X U|| =0), and the same with the fields g and 
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s. From Eq. (.55) one finds that gj 
Euler's equation (.54) that: 



and then from The equation satisfied by F can be easily found from Eq. 



(.58) 



Hill 



u_L(k,i) = — -[u±(k,to) 

a{t) 



dr a(T)s(k, r)] 



to 



(.56) 



Hence, the average transversal component of the pecuUar- 
velocity is damped by the background expansion even in 
the presence of pressure (this is not surprising, since pres- 
sure does not generate vorticity for barotropic fluids, i.e., 
P = Pi.Q))i ^-nd only its fluctuations grow due to the noise. 
The longitudinal component is obtained from the conti- 
nuity equation ( |.53| ), while the gravitational acceleration 
is found from Poisson's equation (.55): 



Ull 



m d5 
k^dt 



g = 



^^'5k 



(.57) 



The elimination of the velocity and gravitational acceler- 
ation fields in favor of the density contrast field by means 
of these expressions yields a closed equation for each mode 
(5(k, t) of the density contrast: 



2H 



d6(k, t) 
dt 



C/fe5(k,t) 



1 



a{t) 



ik • Sj| (k, t) 
(.5 



with 



(.59) 



It is clear from this equation that the noise will induce 
correlations on the density field on top of those already 
present in the initial condition. In the linear regime, this 
superposition is linear and weak noise will only induce 
small Gaussian fluctuations of 5 about the deterministic 
solution. Hence, we can restrict our study to the evolution 
of the average field, (5) , which follows the deterministic 



evolution (i.e., Eq. (.58) dropping the noise). This evolu- 



tion is, however, still different from the dust case because 
of the pressure term. In the rest of the appendix we drop 
the noise source from this equation and understand that 
S really means (5), and the same with the other fields. 

We are particularly interested in the justification of 
the parallelism condition, Eq. (^. From Eqs. {-bT) one 
can easily find a relationship between uy and g: 



g(k, t) ^ F{k,t) U| 



(.60) 



where 



dF 
'dt 



AnGgb ~ HF 



Uk 



■F' 



(.62) 



(i.e., a Ricatti equation). There is parallelism in position 
space between u and g if F does not depend on k; a nec- 
essary condition is that Uk be also fc-independent (i.e., 
there is no pressure) . Therefore, not surprisingly, pressure 
destroys parallelism, but the deviations decrease as one 
goes to larger scales. 

The exact solution to Eq. ( |.62D is difficult to obtain. 
For polytropic models in an Einstein-de Sitter background 
{K = 0, A = 0) one can solve Eq. ( [5^ ) in terms of Meijer's 
G-functions (Haubold et al. 1991) and then compute F 
from its definition, Eq. ( |.61[ ). This procedure is, however, 
algebraically cumbersome and not very illuminating from 
the physical point of view (for some particular values of 
the polytropic index 7 however the solution can be written 
in terms of elementary functions). 

For the purposes of the work presented here, we adopt 
a different approach. We first define Jeans' length Lj by 
0, yielding 



the condition U^-i 



p'iSb) 



-,1/2 



AnGa'^gb 



(.63) 



From Eq. (.58) one can easily grasp the physical mean- 
ing of this quantity: density perturbations with k > Lj^ 
are damped in the linear regime by both pressure and ex- 
pansion, while those with k < Lj^ are damped only by 
expansion (self-gravity dominates over pressure). 

In the limit of small pressure, Lj — » 0, and the pres- 
sure can be considered a perturbation on scales k <C Lj^. 
Defining the small parameter e := L^k'^ (where Li is 
Jeans' length evaluated at some initial time U), we can 
write 



Uk 



AnGgi 



a p'iQb) 
a-i p'{gi) 



(.64) 



and then try a perturbative expansion in powers of e. 
In particular, for polytropic models one has S{t) = 
[a{t) I aiY~^^ , that is, the perturbative expansion is better 
at later times if 7 > 4/3. It should also be noticed that the 
correction due to pressure at a fixed time becomes smaller 
as the polytropic index 7 grows. 

We therefore write -F(k, t) = Yl^=o £"-Pn(k, t) with the 
initial conditions Jo(k, ii) = F(k, ti) and i^„(k, t;) = 0, 
n > 0. From Eq. (|.62|) one finds 



— — = AirGgb 
at 



HFn - F' 



.65) 



F{k,t) = AnGgbit) S{k,t) 



dSjk, t) 
dt 



(.61) 



-dt--^^ 



2Fo)Fi 



sfS 



(.66) 
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Fq obeys the equation of the weU~studied pressureless case 
(Peebles 1980). As i oo it becomes independent of k 
(i.e., it forgets initial conditions) and there is parallelism in 
position space. The asymptotic form of Fq is to be identi- 
fied with F{t) in the parallelism assumption (Eq. [|| in the 
main text), and turns out to be given as: Fq — AnGgbb/b, 



where b(t) is the growing solution of Eq. (.58) when par- 
ticularized to the dust case, p = 0. Fi can be computed 
once Fq is known: 

Fi(k,i) = j'dT S{T)F^{k,T)e- -^'''' [ff(«)+2Fo(k.e)] 

(.67) 

Notice that the whole dependence of Fi on k stems from 
Fq (and hence from initial conditions, which is true for 
every Fn too). The study of this expression can be carried 
out analytically for polytropic models, p — kq^ , in an 
Einstein-de Sitter background. In this case the function 

Fq is 



Fo(k,t) 



t 



5/3 



■ 2A(k) 



i8/3 - 3A(k)t 



(.68) 



with A(k) given by the initial conditions. For polytropic 
models the function S{t) comes down to the simple ex- 
pression S{t) = {t/Uf/^-^f. Taking these results into Eq. 
(.67) then yields 



Fi(k,i) 



-1+27^2/3 



[<5/3 _ 3A(k)]2 



13-67 



¥-27 



67 



1-27 



(.69) 



with the understanding that terms like (t° —t1)/a should 
be substituted by \og{t/ti) when a = 0. 

As t ^ +00, one has the asymptotic behaviors 
Fo(k,i) t-^ and 



Fi(k,t) 



(13-67)i, \t 



t \ 3 



-27 



(7^f), (.70) 



fast enough when k ^ 00, we can Fourier transform this 
expression back to position space and thus obtain the fol- 
lowing relationship between gravitational acceleration and 
the potential component of the velocity. 



g«Fo(i) U|| -L2Fi(i)V2u|| 
which is the result quoted in Sect. |^. 



.73) 
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